3.1.2 \(\int (e x)^m (a+b x^2)^2 (A+B x^2) (c+d x^2) \, dx\) [2]

3.1.2.1 Optimal result
3.1.2.2 Mathematica [A] (verified)
3.1.2.3 Rubi [A] (verified)
3.1.2.4 Maple [B] (verified)
3.1.2.5 Fricas [B] (verification not implemented)
3.1.2.6 Sympy [B] (verification not implemented)
3.1.2.7 Maxima [A] (verification not implemented)
3.1.2.8 Giac [B] (verification not implemented)
3.1.2.9 Mupad [B] (verification not implemented)

3.1.2.1 Optimal result

Integrand size = 29, antiderivative size = 144 \[ \int (e x)^m \left (a+b x^2\right )^2 \left (A+B x^2\right ) \left (c+d x^2\right ) \, dx=\frac {a^2 A c (e x)^{1+m}}{e (1+m)}+\frac {a (2 A b c+a B c+a A d) (e x)^{3+m}}{e^3 (3+m)}+\frac {(a B (2 b c+a d)+A b (b c+2 a d)) (e x)^{5+m}}{e^5 (5+m)}+\frac {b (b B c+A b d+2 a B d) (e x)^{7+m}}{e^7 (7+m)}+\frac {b^2 B d (e x)^{9+m}}{e^9 (9+m)} \]

output
a^2*A*c*(e*x)^(1+m)/e/(1+m)+a*(A*a*d+2*A*b*c+B*a*c)*(e*x)^(3+m)/e^3/(3+m)+ 
(a*B*(a*d+2*b*c)+A*b*(2*a*d+b*c))*(e*x)^(5+m)/e^5/(5+m)+b*(A*b*d+2*B*a*d+B 
*b*c)*(e*x)^(7+m)/e^7/(7+m)+b^2*B*d*(e*x)^(9+m)/e^9/(9+m)
 
3.1.2.2 Mathematica [A] (verified)

Time = 0.13 (sec) , antiderivative size = 113, normalized size of antiderivative = 0.78 \[ \int (e x)^m \left (a+b x^2\right )^2 \left (A+B x^2\right ) \left (c+d x^2\right ) \, dx=x (e x)^m \left (\frac {a^2 A c}{1+m}+\frac {a (2 A b c+a B c+a A d) x^2}{3+m}+\frac {(a B (2 b c+a d)+A b (b c+2 a d)) x^4}{5+m}+\frac {b (b B c+A b d+2 a B d) x^6}{7+m}+\frac {b^2 B d x^8}{9+m}\right ) \]

input
Integrate[(e*x)^m*(a + b*x^2)^2*(A + B*x^2)*(c + d*x^2),x]
 
output
x*(e*x)^m*((a^2*A*c)/(1 + m) + (a*(2*A*b*c + a*B*c + a*A*d)*x^2)/(3 + m) + 
 ((a*B*(2*b*c + a*d) + A*b*(b*c + 2*a*d))*x^4)/(5 + m) + (b*(b*B*c + A*b*d 
 + 2*a*B*d)*x^6)/(7 + m) + (b^2*B*d*x^8)/(9 + m))
 
3.1.2.3 Rubi [A] (verified)

Time = 0.32 (sec) , antiderivative size = 144, normalized size of antiderivative = 1.00, number of steps used = 2, number of rules used = 2, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.069, Rules used = {437, 2009}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int \left (a+b x^2\right )^2 \left (A+B x^2\right ) \left (c+d x^2\right ) (e x)^m \, dx\)

\(\Big \downarrow \) 437

\(\displaystyle \int \left (a^2 A c (e x)^m+\frac {b (e x)^{m+6} (2 a B d+A b d+b B c)}{e^6}+\frac {(e x)^{m+4} (A b (2 a d+b c)+a B (a d+2 b c))}{e^4}+\frac {a (e x)^{m+2} (a A d+a B c+2 A b c)}{e^2}+\frac {b^2 B d (e x)^{m+8}}{e^8}\right )dx\)

\(\Big \downarrow \) 2009

\(\displaystyle \frac {a^2 A c (e x)^{m+1}}{e (m+1)}+\frac {b (e x)^{m+7} (2 a B d+A b d+b B c)}{e^7 (m+7)}+\frac {(e x)^{m+5} (A b (2 a d+b c)+a B (a d+2 b c))}{e^5 (m+5)}+\frac {a (e x)^{m+3} (a A d+a B c+2 A b c)}{e^3 (m+3)}+\frac {b^2 B d (e x)^{m+9}}{e^9 (m+9)}\)

input
Int[(e*x)^m*(a + b*x^2)^2*(A + B*x^2)*(c + d*x^2),x]
 
output
(a^2*A*c*(e*x)^(1 + m))/(e*(1 + m)) + (a*(2*A*b*c + a*B*c + a*A*d)*(e*x)^( 
3 + m))/(e^3*(3 + m)) + ((a*B*(2*b*c + a*d) + A*b*(b*c + 2*a*d))*(e*x)^(5 
+ m))/(e^5*(5 + m)) + (b*(b*B*c + A*b*d + 2*a*B*d)*(e*x)^(7 + m))/(e^7*(7 
+ m)) + (b^2*B*d*(e*x)^(9 + m))/(e^9*(9 + m))
 

3.1.2.3.1 Defintions of rubi rules used

rule 437
Int[((g_.)*(x_))^(m_.)*((a_) + (b_.)*(x_)^2)^(p_.)*((c_) + (d_.)*(x_)^2)^(q 
_.)*((e_) + (f_.)*(x_)^2)^(r_.), x_Symbol] :> Int[ExpandIntegrand[(g*x)^m*( 
a + b*x^2)^p*(c + d*x^2)^q*(e + f*x^2)^r, x], x] /; FreeQ[{a, b, c, d, e, f 
, g, m}, x] && IGtQ[p, -2] && IGtQ[q, 0] && IGtQ[r, 0]
 

rule 2009
Int[u_, x_Symbol] :> Simp[IntSum[u, x], x] /; SumQ[u]
 
3.1.2.4 Maple [B] (verified)

Leaf count of result is larger than twice the leaf count of optimal. \(710\) vs. \(2(144)=288\).

Time = 3.38 (sec) , antiderivative size = 711, normalized size of antiderivative = 4.94

method result size
gosper \(\frac {x \left (B \,b^{2} d \,m^{4} x^{8}+16 B \,b^{2} d \,m^{3} x^{8}+A \,b^{2} d \,m^{4} x^{6}+2 B a b d \,m^{4} x^{6}+B \,b^{2} c \,m^{4} x^{6}+86 B \,b^{2} d \,m^{2} x^{8}+18 A \,b^{2} d \,m^{3} x^{6}+36 B a b d \,m^{3} x^{6}+18 B \,b^{2} c \,m^{3} x^{6}+176 m \,x^{8} b^{2} B d +2 A a b d \,m^{4} x^{4}+A \,b^{2} c \,m^{4} x^{4}+104 A \,b^{2} d \,m^{2} x^{6}+B \,a^{2} d \,m^{4} x^{4}+2 B a b c \,m^{4} x^{4}+208 B a b d \,m^{2} x^{6}+104 B \,b^{2} c \,m^{2} x^{6}+105 b^{2} B d \,x^{8}+40 A a b d \,m^{3} x^{4}+20 A \,b^{2} c \,m^{3} x^{4}+222 A \,b^{2} d \,x^{6} m +20 B \,a^{2} d \,m^{3} x^{4}+40 B a b c \,m^{3} x^{4}+444 B a b d \,x^{6} m +222 B \,b^{2} c \,x^{6} m +A \,a^{2} d \,m^{4} x^{2}+2 A a b c \,m^{4} x^{2}+260 A a b d \,m^{2} x^{4}+130 A \,b^{2} c \,m^{2} x^{4}+135 A \,b^{2} d \,x^{6}+B \,a^{2} c \,m^{4} x^{2}+130 B \,a^{2} d \,m^{2} x^{4}+260 B a b c \,m^{2} x^{4}+270 B a b d \,x^{6}+135 B \,b^{2} c \,x^{6}+22 A \,a^{2} d \,m^{3} x^{2}+44 A a b c \,m^{3} x^{2}+600 A a b d \,x^{4} m +300 A \,b^{2} c \,x^{4} m +22 B \,a^{2} c \,m^{3} x^{2}+300 B \,a^{2} d \,x^{4} m +600 B a b c \,x^{4} m +A \,a^{2} c \,m^{4}+164 A \,a^{2} d \,m^{2} x^{2}+328 A a b c \,m^{2} x^{2}+378 A a b d \,x^{4}+189 A \,b^{2} c \,x^{4}+164 B \,a^{2} c \,m^{2} x^{2}+189 B \,a^{2} d \,x^{4}+378 B a b c \,x^{4}+24 A \,a^{2} c \,m^{3}+458 A \,a^{2} d \,x^{2} m +916 A a b c \,x^{2} m +458 B \,a^{2} c \,x^{2} m +206 A \,a^{2} c \,m^{2}+315 A \,a^{2} d \,x^{2}+630 A a b c \,x^{2}+315 B \,a^{2} c \,x^{2}+744 A \,a^{2} c m +945 A \,a^{2} c \right ) \left (e x \right )^{m}}{\left (9+m \right ) \left (7+m \right ) \left (5+m \right ) \left (3+m \right ) \left (1+m \right )}\) \(711\)
risch \(\frac {x \left (B \,b^{2} d \,m^{4} x^{8}+16 B \,b^{2} d \,m^{3} x^{8}+A \,b^{2} d \,m^{4} x^{6}+2 B a b d \,m^{4} x^{6}+B \,b^{2} c \,m^{4} x^{6}+86 B \,b^{2} d \,m^{2} x^{8}+18 A \,b^{2} d \,m^{3} x^{6}+36 B a b d \,m^{3} x^{6}+18 B \,b^{2} c \,m^{3} x^{6}+176 m \,x^{8} b^{2} B d +2 A a b d \,m^{4} x^{4}+A \,b^{2} c \,m^{4} x^{4}+104 A \,b^{2} d \,m^{2} x^{6}+B \,a^{2} d \,m^{4} x^{4}+2 B a b c \,m^{4} x^{4}+208 B a b d \,m^{2} x^{6}+104 B \,b^{2} c \,m^{2} x^{6}+105 b^{2} B d \,x^{8}+40 A a b d \,m^{3} x^{4}+20 A \,b^{2} c \,m^{3} x^{4}+222 A \,b^{2} d \,x^{6} m +20 B \,a^{2} d \,m^{3} x^{4}+40 B a b c \,m^{3} x^{4}+444 B a b d \,x^{6} m +222 B \,b^{2} c \,x^{6} m +A \,a^{2} d \,m^{4} x^{2}+2 A a b c \,m^{4} x^{2}+260 A a b d \,m^{2} x^{4}+130 A \,b^{2} c \,m^{2} x^{4}+135 A \,b^{2} d \,x^{6}+B \,a^{2} c \,m^{4} x^{2}+130 B \,a^{2} d \,m^{2} x^{4}+260 B a b c \,m^{2} x^{4}+270 B a b d \,x^{6}+135 B \,b^{2} c \,x^{6}+22 A \,a^{2} d \,m^{3} x^{2}+44 A a b c \,m^{3} x^{2}+600 A a b d \,x^{4} m +300 A \,b^{2} c \,x^{4} m +22 B \,a^{2} c \,m^{3} x^{2}+300 B \,a^{2} d \,x^{4} m +600 B a b c \,x^{4} m +A \,a^{2} c \,m^{4}+164 A \,a^{2} d \,m^{2} x^{2}+328 A a b c \,m^{2} x^{2}+378 A a b d \,x^{4}+189 A \,b^{2} c \,x^{4}+164 B \,a^{2} c \,m^{2} x^{2}+189 B \,a^{2} d \,x^{4}+378 B a b c \,x^{4}+24 A \,a^{2} c \,m^{3}+458 A \,a^{2} d \,x^{2} m +916 A a b c \,x^{2} m +458 B \,a^{2} c \,x^{2} m +206 A \,a^{2} c \,m^{2}+315 A \,a^{2} d \,x^{2}+630 A a b c \,x^{2}+315 B \,a^{2} c \,x^{2}+744 A \,a^{2} c m +945 A \,a^{2} c \right ) \left (e x \right )^{m}}{\left (9+m \right ) \left (7+m \right ) \left (5+m \right ) \left (3+m \right ) \left (1+m \right )}\) \(711\)
parallelrisch \(\text {Expression too large to display}\) \(1010\)

input
int((e*x)^m*(b*x^2+a)^2*(B*x^2+A)*(d*x^2+c),x,method=_RETURNVERBOSE)
 
output
x*(B*b^2*d*m^4*x^8+16*B*b^2*d*m^3*x^8+A*b^2*d*m^4*x^6+2*B*a*b*d*m^4*x^6+B* 
b^2*c*m^4*x^6+86*B*b^2*d*m^2*x^8+18*A*b^2*d*m^3*x^6+36*B*a*b*d*m^3*x^6+18* 
B*b^2*c*m^3*x^6+176*B*b^2*d*m*x^8+2*A*a*b*d*m^4*x^4+A*b^2*c*m^4*x^4+104*A* 
b^2*d*m^2*x^6+B*a^2*d*m^4*x^4+2*B*a*b*c*m^4*x^4+208*B*a*b*d*m^2*x^6+104*B* 
b^2*c*m^2*x^6+105*B*b^2*d*x^8+40*A*a*b*d*m^3*x^4+20*A*b^2*c*m^3*x^4+222*A* 
b^2*d*m*x^6+20*B*a^2*d*m^3*x^4+40*B*a*b*c*m^3*x^4+444*B*a*b*d*m*x^6+222*B* 
b^2*c*m*x^6+A*a^2*d*m^4*x^2+2*A*a*b*c*m^4*x^2+260*A*a*b*d*m^2*x^4+130*A*b^ 
2*c*m^2*x^4+135*A*b^2*d*x^6+B*a^2*c*m^4*x^2+130*B*a^2*d*m^2*x^4+260*B*a*b* 
c*m^2*x^4+270*B*a*b*d*x^6+135*B*b^2*c*x^6+22*A*a^2*d*m^3*x^2+44*A*a*b*c*m^ 
3*x^2+600*A*a*b*d*m*x^4+300*A*b^2*c*m*x^4+22*B*a^2*c*m^3*x^2+300*B*a^2*d*m 
*x^4+600*B*a*b*c*m*x^4+A*a^2*c*m^4+164*A*a^2*d*m^2*x^2+328*A*a*b*c*m^2*x^2 
+378*A*a*b*d*x^4+189*A*b^2*c*x^4+164*B*a^2*c*m^2*x^2+189*B*a^2*d*x^4+378*B 
*a*b*c*x^4+24*A*a^2*c*m^3+458*A*a^2*d*m*x^2+916*A*a*b*c*m*x^2+458*B*a^2*c* 
m*x^2+206*A*a^2*c*m^2+315*A*a^2*d*x^2+630*A*a*b*c*x^2+315*B*a^2*c*x^2+744* 
A*a^2*c*m+945*A*a^2*c)*(e*x)^m/(9+m)/(7+m)/(5+m)/(3+m)/(1+m)
 
3.1.2.5 Fricas [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 532 vs. \(2 (144) = 288\).

Time = 0.28 (sec) , antiderivative size = 532, normalized size of antiderivative = 3.69 \[ \int (e x)^m \left (a+b x^2\right )^2 \left (A+B x^2\right ) \left (c+d x^2\right ) \, dx=\frac {{\left ({\left (B b^{2} d m^{4} + 16 \, B b^{2} d m^{3} + 86 \, B b^{2} d m^{2} + 176 \, B b^{2} d m + 105 \, B b^{2} d\right )} x^{9} + {\left ({\left (B b^{2} c + {\left (2 \, B a b + A b^{2}\right )} d\right )} m^{4} + 135 \, B b^{2} c + 18 \, {\left (B b^{2} c + {\left (2 \, B a b + A b^{2}\right )} d\right )} m^{3} + 104 \, {\left (B b^{2} c + {\left (2 \, B a b + A b^{2}\right )} d\right )} m^{2} + 135 \, {\left (2 \, B a b + A b^{2}\right )} d + 222 \, {\left (B b^{2} c + {\left (2 \, B a b + A b^{2}\right )} d\right )} m\right )} x^{7} + {\left ({\left ({\left (2 \, B a b + A b^{2}\right )} c + {\left (B a^{2} + 2 \, A a b\right )} d\right )} m^{4} + 20 \, {\left ({\left (2 \, B a b + A b^{2}\right )} c + {\left (B a^{2} + 2 \, A a b\right )} d\right )} m^{3} + 130 \, {\left ({\left (2 \, B a b + A b^{2}\right )} c + {\left (B a^{2} + 2 \, A a b\right )} d\right )} m^{2} + 189 \, {\left (2 \, B a b + A b^{2}\right )} c + 189 \, {\left (B a^{2} + 2 \, A a b\right )} d + 300 \, {\left ({\left (2 \, B a b + A b^{2}\right )} c + {\left (B a^{2} + 2 \, A a b\right )} d\right )} m\right )} x^{5} + {\left ({\left (A a^{2} d + {\left (B a^{2} + 2 \, A a b\right )} c\right )} m^{4} + 315 \, A a^{2} d + 22 \, {\left (A a^{2} d + {\left (B a^{2} + 2 \, A a b\right )} c\right )} m^{3} + 164 \, {\left (A a^{2} d + {\left (B a^{2} + 2 \, A a b\right )} c\right )} m^{2} + 315 \, {\left (B a^{2} + 2 \, A a b\right )} c + 458 \, {\left (A a^{2} d + {\left (B a^{2} + 2 \, A a b\right )} c\right )} m\right )} x^{3} + {\left (A a^{2} c m^{4} + 24 \, A a^{2} c m^{3} + 206 \, A a^{2} c m^{2} + 744 \, A a^{2} c m + 945 \, A a^{2} c\right )} x\right )} \left (e x\right )^{m}}{m^{5} + 25 \, m^{4} + 230 \, m^{3} + 950 \, m^{2} + 1689 \, m + 945} \]

input
integrate((e*x)^m*(b*x^2+a)^2*(B*x^2+A)*(d*x^2+c),x, algorithm="fricas")
 
output
((B*b^2*d*m^4 + 16*B*b^2*d*m^3 + 86*B*b^2*d*m^2 + 176*B*b^2*d*m + 105*B*b^ 
2*d)*x^9 + ((B*b^2*c + (2*B*a*b + A*b^2)*d)*m^4 + 135*B*b^2*c + 18*(B*b^2* 
c + (2*B*a*b + A*b^2)*d)*m^3 + 104*(B*b^2*c + (2*B*a*b + A*b^2)*d)*m^2 + 1 
35*(2*B*a*b + A*b^2)*d + 222*(B*b^2*c + (2*B*a*b + A*b^2)*d)*m)*x^7 + (((2 
*B*a*b + A*b^2)*c + (B*a^2 + 2*A*a*b)*d)*m^4 + 20*((2*B*a*b + A*b^2)*c + ( 
B*a^2 + 2*A*a*b)*d)*m^3 + 130*((2*B*a*b + A*b^2)*c + (B*a^2 + 2*A*a*b)*d)* 
m^2 + 189*(2*B*a*b + A*b^2)*c + 189*(B*a^2 + 2*A*a*b)*d + 300*((2*B*a*b + 
A*b^2)*c + (B*a^2 + 2*A*a*b)*d)*m)*x^5 + ((A*a^2*d + (B*a^2 + 2*A*a*b)*c)* 
m^4 + 315*A*a^2*d + 22*(A*a^2*d + (B*a^2 + 2*A*a*b)*c)*m^3 + 164*(A*a^2*d 
+ (B*a^2 + 2*A*a*b)*c)*m^2 + 315*(B*a^2 + 2*A*a*b)*c + 458*(A*a^2*d + (B*a 
^2 + 2*A*a*b)*c)*m)*x^3 + (A*a^2*c*m^4 + 24*A*a^2*c*m^3 + 206*A*a^2*c*m^2 
+ 744*A*a^2*c*m + 945*A*a^2*c)*x)*(e*x)^m/(m^5 + 25*m^4 + 230*m^3 + 950*m^ 
2 + 1689*m + 945)
 
3.1.2.6 Sympy [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 3271 vs. \(2 (139) = 278\).

Time = 0.66 (sec) , antiderivative size = 3271, normalized size of antiderivative = 22.72 \[ \int (e x)^m \left (a+b x^2\right )^2 \left (A+B x^2\right ) \left (c+d x^2\right ) \, dx=\text {Too large to display} \]

input
integrate((e*x)**m*(b*x**2+a)**2*(B*x**2+A)*(d*x**2+c),x)
 
output
Piecewise(((-A*a**2*c/(8*x**8) - A*a**2*d/(6*x**6) - A*a*b*c/(3*x**6) - A* 
a*b*d/(2*x**4) - A*b**2*c/(4*x**4) - A*b**2*d/(2*x**2) - B*a**2*c/(6*x**6) 
 - B*a**2*d/(4*x**4) - B*a*b*c/(2*x**4) - B*a*b*d/x**2 - B*b**2*c/(2*x**2) 
 + B*b**2*d*log(x))/e**9, Eq(m, -9)), ((-A*a**2*c/(6*x**6) - A*a**2*d/(4*x 
**4) - A*a*b*c/(2*x**4) - A*a*b*d/x**2 - A*b**2*c/(2*x**2) + A*b**2*d*log( 
x) - B*a**2*c/(4*x**4) - B*a**2*d/(2*x**2) - B*a*b*c/x**2 + 2*B*a*b*d*log( 
x) + B*b**2*c*log(x) + B*b**2*d*x**2/2)/e**7, Eq(m, -7)), ((-A*a**2*c/(4*x 
**4) - A*a**2*d/(2*x**2) - A*a*b*c/x**2 + 2*A*a*b*d*log(x) + A*b**2*c*log( 
x) + A*b**2*d*x**2/2 - B*a**2*c/(2*x**2) + B*a**2*d*log(x) + 2*B*a*b*c*log 
(x) + B*a*b*d*x**2 + B*b**2*c*x**2/2 + B*b**2*d*x**4/4)/e**5, Eq(m, -5)), 
((-A*a**2*c/(2*x**2) + A*a**2*d*log(x) + 2*A*a*b*c*log(x) + A*a*b*d*x**2 + 
 A*b**2*c*x**2/2 + A*b**2*d*x**4/4 + B*a**2*c*log(x) + B*a**2*d*x**2/2 + B 
*a*b*c*x**2 + B*a*b*d*x**4/2 + B*b**2*c*x**4/4 + B*b**2*d*x**6/6)/e**3, Eq 
(m, -3)), ((A*a**2*c*log(x) + A*a**2*d*x**2/2 + A*a*b*c*x**2 + A*a*b*d*x** 
4/2 + A*b**2*c*x**4/4 + A*b**2*d*x**6/6 + B*a**2*c*x**2/2 + B*a**2*d*x**4/ 
4 + B*a*b*c*x**4/2 + B*a*b*d*x**6/3 + B*b**2*c*x**6/6 + B*b**2*d*x**8/8)/e 
, Eq(m, -1)), (A*a**2*c*m**4*x*(e*x)**m/(m**5 + 25*m**4 + 230*m**3 + 950*m 
**2 + 1689*m + 945) + 24*A*a**2*c*m**3*x*(e*x)**m/(m**5 + 25*m**4 + 230*m* 
*3 + 950*m**2 + 1689*m + 945) + 206*A*a**2*c*m**2*x*(e*x)**m/(m**5 + 25*m* 
*4 + 230*m**3 + 950*m**2 + 1689*m + 945) + 744*A*a**2*c*m*x*(e*x)**m/(m...
 
3.1.2.7 Maxima [A] (verification not implemented)

Time = 0.23 (sec) , antiderivative size = 242, normalized size of antiderivative = 1.68 \[ \int (e x)^m \left (a+b x^2\right )^2 \left (A+B x^2\right ) \left (c+d x^2\right ) \, dx=\frac {B b^{2} d e^{m} x^{9} x^{m}}{m + 9} + \frac {B b^{2} c e^{m} x^{7} x^{m}}{m + 7} + \frac {2 \, B a b d e^{m} x^{7} x^{m}}{m + 7} + \frac {A b^{2} d e^{m} x^{7} x^{m}}{m + 7} + \frac {2 \, B a b c e^{m} x^{5} x^{m}}{m + 5} + \frac {A b^{2} c e^{m} x^{5} x^{m}}{m + 5} + \frac {B a^{2} d e^{m} x^{5} x^{m}}{m + 5} + \frac {2 \, A a b d e^{m} x^{5} x^{m}}{m + 5} + \frac {B a^{2} c e^{m} x^{3} x^{m}}{m + 3} + \frac {2 \, A a b c e^{m} x^{3} x^{m}}{m + 3} + \frac {A a^{2} d e^{m} x^{3} x^{m}}{m + 3} + \frac {\left (e x\right )^{m + 1} A a^{2} c}{e {\left (m + 1\right )}} \]

input
integrate((e*x)^m*(b*x^2+a)^2*(B*x^2+A)*(d*x^2+c),x, algorithm="maxima")
 
output
B*b^2*d*e^m*x^9*x^m/(m + 9) + B*b^2*c*e^m*x^7*x^m/(m + 7) + 2*B*a*b*d*e^m* 
x^7*x^m/(m + 7) + A*b^2*d*e^m*x^7*x^m/(m + 7) + 2*B*a*b*c*e^m*x^5*x^m/(m + 
 5) + A*b^2*c*e^m*x^5*x^m/(m + 5) + B*a^2*d*e^m*x^5*x^m/(m + 5) + 2*A*a*b* 
d*e^m*x^5*x^m/(m + 5) + B*a^2*c*e^m*x^3*x^m/(m + 3) + 2*A*a*b*c*e^m*x^3*x^ 
m/(m + 3) + A*a^2*d*e^m*x^3*x^m/(m + 3) + (e*x)^(m + 1)*A*a^2*c/(e*(m + 1) 
)
 
3.1.2.8 Giac [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 1009 vs. \(2 (144) = 288\).

Time = 0.30 (sec) , antiderivative size = 1009, normalized size of antiderivative = 7.01 \[ \int (e x)^m \left (a+b x^2\right )^2 \left (A+B x^2\right ) \left (c+d x^2\right ) \, dx=\frac {\left (e x\right )^{m} B b^{2} d m^{4} x^{9} + 16 \, \left (e x\right )^{m} B b^{2} d m^{3} x^{9} + \left (e x\right )^{m} B b^{2} c m^{4} x^{7} + 2 \, \left (e x\right )^{m} B a b d m^{4} x^{7} + \left (e x\right )^{m} A b^{2} d m^{4} x^{7} + 86 \, \left (e x\right )^{m} B b^{2} d m^{2} x^{9} + 18 \, \left (e x\right )^{m} B b^{2} c m^{3} x^{7} + 36 \, \left (e x\right )^{m} B a b d m^{3} x^{7} + 18 \, \left (e x\right )^{m} A b^{2} d m^{3} x^{7} + 176 \, \left (e x\right )^{m} B b^{2} d m x^{9} + 2 \, \left (e x\right )^{m} B a b c m^{4} x^{5} + \left (e x\right )^{m} A b^{2} c m^{4} x^{5} + \left (e x\right )^{m} B a^{2} d m^{4} x^{5} + 2 \, \left (e x\right )^{m} A a b d m^{4} x^{5} + 104 \, \left (e x\right )^{m} B b^{2} c m^{2} x^{7} + 208 \, \left (e x\right )^{m} B a b d m^{2} x^{7} + 104 \, \left (e x\right )^{m} A b^{2} d m^{2} x^{7} + 105 \, \left (e x\right )^{m} B b^{2} d x^{9} + 40 \, \left (e x\right )^{m} B a b c m^{3} x^{5} + 20 \, \left (e x\right )^{m} A b^{2} c m^{3} x^{5} + 20 \, \left (e x\right )^{m} B a^{2} d m^{3} x^{5} + 40 \, \left (e x\right )^{m} A a b d m^{3} x^{5} + 222 \, \left (e x\right )^{m} B b^{2} c m x^{7} + 444 \, \left (e x\right )^{m} B a b d m x^{7} + 222 \, \left (e x\right )^{m} A b^{2} d m x^{7} + \left (e x\right )^{m} B a^{2} c m^{4} x^{3} + 2 \, \left (e x\right )^{m} A a b c m^{4} x^{3} + \left (e x\right )^{m} A a^{2} d m^{4} x^{3} + 260 \, \left (e x\right )^{m} B a b c m^{2} x^{5} + 130 \, \left (e x\right )^{m} A b^{2} c m^{2} x^{5} + 130 \, \left (e x\right )^{m} B a^{2} d m^{2} x^{5} + 260 \, \left (e x\right )^{m} A a b d m^{2} x^{5} + 135 \, \left (e x\right )^{m} B b^{2} c x^{7} + 270 \, \left (e x\right )^{m} B a b d x^{7} + 135 \, \left (e x\right )^{m} A b^{2} d x^{7} + 22 \, \left (e x\right )^{m} B a^{2} c m^{3} x^{3} + 44 \, \left (e x\right )^{m} A a b c m^{3} x^{3} + 22 \, \left (e x\right )^{m} A a^{2} d m^{3} x^{3} + 600 \, \left (e x\right )^{m} B a b c m x^{5} + 300 \, \left (e x\right )^{m} A b^{2} c m x^{5} + 300 \, \left (e x\right )^{m} B a^{2} d m x^{5} + 600 \, \left (e x\right )^{m} A a b d m x^{5} + \left (e x\right )^{m} A a^{2} c m^{4} x + 164 \, \left (e x\right )^{m} B a^{2} c m^{2} x^{3} + 328 \, \left (e x\right )^{m} A a b c m^{2} x^{3} + 164 \, \left (e x\right )^{m} A a^{2} d m^{2} x^{3} + 378 \, \left (e x\right )^{m} B a b c x^{5} + 189 \, \left (e x\right )^{m} A b^{2} c x^{5} + 189 \, \left (e x\right )^{m} B a^{2} d x^{5} + 378 \, \left (e x\right )^{m} A a b d x^{5} + 24 \, \left (e x\right )^{m} A a^{2} c m^{3} x + 458 \, \left (e x\right )^{m} B a^{2} c m x^{3} + 916 \, \left (e x\right )^{m} A a b c m x^{3} + 458 \, \left (e x\right )^{m} A a^{2} d m x^{3} + 206 \, \left (e x\right )^{m} A a^{2} c m^{2} x + 315 \, \left (e x\right )^{m} B a^{2} c x^{3} + 630 \, \left (e x\right )^{m} A a b c x^{3} + 315 \, \left (e x\right )^{m} A a^{2} d x^{3} + 744 \, \left (e x\right )^{m} A a^{2} c m x + 945 \, \left (e x\right )^{m} A a^{2} c x}{m^{5} + 25 \, m^{4} + 230 \, m^{3} + 950 \, m^{2} + 1689 \, m + 945} \]

input
integrate((e*x)^m*(b*x^2+a)^2*(B*x^2+A)*(d*x^2+c),x, algorithm="giac")
 
output
((e*x)^m*B*b^2*d*m^4*x^9 + 16*(e*x)^m*B*b^2*d*m^3*x^9 + (e*x)^m*B*b^2*c*m^ 
4*x^7 + 2*(e*x)^m*B*a*b*d*m^4*x^7 + (e*x)^m*A*b^2*d*m^4*x^7 + 86*(e*x)^m*B 
*b^2*d*m^2*x^9 + 18*(e*x)^m*B*b^2*c*m^3*x^7 + 36*(e*x)^m*B*a*b*d*m^3*x^7 + 
 18*(e*x)^m*A*b^2*d*m^3*x^7 + 176*(e*x)^m*B*b^2*d*m*x^9 + 2*(e*x)^m*B*a*b* 
c*m^4*x^5 + (e*x)^m*A*b^2*c*m^4*x^5 + (e*x)^m*B*a^2*d*m^4*x^5 + 2*(e*x)^m* 
A*a*b*d*m^4*x^5 + 104*(e*x)^m*B*b^2*c*m^2*x^7 + 208*(e*x)^m*B*a*b*d*m^2*x^ 
7 + 104*(e*x)^m*A*b^2*d*m^2*x^7 + 105*(e*x)^m*B*b^2*d*x^9 + 40*(e*x)^m*B*a 
*b*c*m^3*x^5 + 20*(e*x)^m*A*b^2*c*m^3*x^5 + 20*(e*x)^m*B*a^2*d*m^3*x^5 + 4 
0*(e*x)^m*A*a*b*d*m^3*x^5 + 222*(e*x)^m*B*b^2*c*m*x^7 + 444*(e*x)^m*B*a*b* 
d*m*x^7 + 222*(e*x)^m*A*b^2*d*m*x^7 + (e*x)^m*B*a^2*c*m^4*x^3 + 2*(e*x)^m* 
A*a*b*c*m^4*x^3 + (e*x)^m*A*a^2*d*m^4*x^3 + 260*(e*x)^m*B*a*b*c*m^2*x^5 + 
130*(e*x)^m*A*b^2*c*m^2*x^5 + 130*(e*x)^m*B*a^2*d*m^2*x^5 + 260*(e*x)^m*A* 
a*b*d*m^2*x^5 + 135*(e*x)^m*B*b^2*c*x^7 + 270*(e*x)^m*B*a*b*d*x^7 + 135*(e 
*x)^m*A*b^2*d*x^7 + 22*(e*x)^m*B*a^2*c*m^3*x^3 + 44*(e*x)^m*A*a*b*c*m^3*x^ 
3 + 22*(e*x)^m*A*a^2*d*m^3*x^3 + 600*(e*x)^m*B*a*b*c*m*x^5 + 300*(e*x)^m*A 
*b^2*c*m*x^5 + 300*(e*x)^m*B*a^2*d*m*x^5 + 600*(e*x)^m*A*a*b*d*m*x^5 + (e* 
x)^m*A*a^2*c*m^4*x + 164*(e*x)^m*B*a^2*c*m^2*x^3 + 328*(e*x)^m*A*a*b*c*m^2 
*x^3 + 164*(e*x)^m*A*a^2*d*m^2*x^3 + 378*(e*x)^m*B*a*b*c*x^5 + 189*(e*x)^m 
*A*b^2*c*x^5 + 189*(e*x)^m*B*a^2*d*x^5 + 378*(e*x)^m*A*a*b*d*x^5 + 24*(e*x 
)^m*A*a^2*c*m^3*x + 458*(e*x)^m*B*a^2*c*m*x^3 + 916*(e*x)^m*A*a*b*c*m*x...
 
3.1.2.9 Mupad [B] (verification not implemented)

Time = 5.63 (sec) , antiderivative size = 305, normalized size of antiderivative = 2.12 \[ \int (e x)^m \left (a+b x^2\right )^2 \left (A+B x^2\right ) \left (c+d x^2\right ) \, dx={\left (e\,x\right )}^m\,\left (\frac {x^5\,\left (A\,b^2\,c+B\,a^2\,d+2\,A\,a\,b\,d+2\,B\,a\,b\,c\right )\,\left (m^4+20\,m^3+130\,m^2+300\,m+189\right )}{m^5+25\,m^4+230\,m^3+950\,m^2+1689\,m+945}+\frac {a\,x^3\,\left (A\,a\,d+2\,A\,b\,c+B\,a\,c\right )\,\left (m^4+22\,m^3+164\,m^2+458\,m+315\right )}{m^5+25\,m^4+230\,m^3+950\,m^2+1689\,m+945}+\frac {b\,x^7\,\left (A\,b\,d+2\,B\,a\,d+B\,b\,c\right )\,\left (m^4+18\,m^3+104\,m^2+222\,m+135\right )}{m^5+25\,m^4+230\,m^3+950\,m^2+1689\,m+945}+\frac {A\,a^2\,c\,x\,\left (m^4+24\,m^3+206\,m^2+744\,m+945\right )}{m^5+25\,m^4+230\,m^3+950\,m^2+1689\,m+945}+\frac {B\,b^2\,d\,x^9\,\left (m^4+16\,m^3+86\,m^2+176\,m+105\right )}{m^5+25\,m^4+230\,m^3+950\,m^2+1689\,m+945}\right ) \]

input
int((A + B*x^2)*(e*x)^m*(a + b*x^2)^2*(c + d*x^2),x)
 
output
(e*x)^m*((x^5*(A*b^2*c + B*a^2*d + 2*A*a*b*d + 2*B*a*b*c)*(300*m + 130*m^2 
 + 20*m^3 + m^4 + 189))/(1689*m + 950*m^2 + 230*m^3 + 25*m^4 + m^5 + 945) 
+ (a*x^3*(A*a*d + 2*A*b*c + B*a*c)*(458*m + 164*m^2 + 22*m^3 + m^4 + 315)) 
/(1689*m + 950*m^2 + 230*m^3 + 25*m^4 + m^5 + 945) + (b*x^7*(A*b*d + 2*B*a 
*d + B*b*c)*(222*m + 104*m^2 + 18*m^3 + m^4 + 135))/(1689*m + 950*m^2 + 23 
0*m^3 + 25*m^4 + m^5 + 945) + (A*a^2*c*x*(744*m + 206*m^2 + 24*m^3 + m^4 + 
 945))/(1689*m + 950*m^2 + 230*m^3 + 25*m^4 + m^5 + 945) + (B*b^2*d*x^9*(1 
76*m + 86*m^2 + 16*m^3 + m^4 + 105))/(1689*m + 950*m^2 + 230*m^3 + 25*m^4 
+ m^5 + 945))